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Abstract. 

We report on a systematic approach for the calculation of the negativity in the 
ground state of a one-dimensional quantum field theory. The partial transpose p^ 2 
of the reduced density matrix of a subsystem A = A\ VJA2 is explicitly constructed 
as an imaginary-time path integral and from this the replicated traces Tr(p^ 2 ) n 
are obtained. The logarithmic negativity S = log | |p^ 2 1 1 is then the continuation to 
n — > 1 of the traces of the even powers. For pure states, this procedure reproduces 
the known results. We then apply this method to conformally invariant field 
theories in several different physical situations for infinite and finite systems and 
without or with boundaries. In particular, in the case of two adjacent intervals of 
lengths £1, ^2 in an infinite system, wc derive the result £ ~ (c/4) ln(^i<?2 / (^1+^2)), 
where c is the central charge. For the more complicated case of two disjoint 
intervals, we show that the negativity depends only on the harmonic ratio of the 
four end-points and so is manifestly scale invariant. Wc explicitly calculate the 
scale-invariant functions for the replicated traces in the case of the CFT for the 
free compactificd boson, but we have not so far been able to obtain the n — » 1 
continuation for the negativity even in the limit of large compactification radius. 
We have checked all our findings against exact numerical results for the harmonic 
chain which is described by a non-compactified free boson. 
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1. Introduction 

The study of the entanglement content of many-body quantum systems has allowed in 
recent years a deeper understanding of these systems, in particular in connection with 
criticality and topological order (see [1] for reviews) . In the case when a system is in its 
ground-state (or any given pure state) the entanglement between two complementary 
parts is measured by the entanglement entropies defined as follows. Let p be the 
density matrix of a system, which we take to be in a pure quantum state \^S), so that 
p = Let the Hilbert space be written as a direct product % = Ha ®Hb- 

reduced density matrix is pa — Tr^p. The entanglement entropy is the corresponding 
von Neumann entropy 

S A = -Tvp A lnp A , (1) 

and analogously for Sb • When p corresponds to a pure quantum state Sa — Sb- Other 
standard measures of bipartite entanglement in pure states are the Renyi entropies 

5 i n) = nb; lnTr ^, (2) 

that also satisfy = S^ whenever p corresponds to a pure quantum state. 

From these definitions Sa = lirn . All for any n and for pure states are 

entanglement monotones [2], i.e. are quantities which do not increase under LOCC 
(local operation and classical communication), which is a key property of any quantity 
to be a good measure of entanglement. The knowledge of the Renyi entropies for any 
n determines also the full spectrum of the reduced density matrix [3]. 

For a mixed state the entanglement entropies are not longer good measures of 
entanglement since they mix quantum and classical correlations (e.g. in an high 
temperature mixed state, Sa gives the extensive result for the thermal entropy that 
has nothing to do with entanglement). This is also evident from the fact that Sa is 
no longer equal to Sb- A quantity that is easily constructed from S^ and Sg^ is the 
(Renyi) mutual information, defined as 

r(«) _ c(") , _ c(") (o\ 

1 A:B~°A ~r D B D AuB > \°) 

that by definition is symmetric in A and B. However, I A n B has not all the correct 
properties to be an entanglement measure and indeed it is not an entanglement 
monotone for any n (for example it has been shown that for some separable mixed 
states it is non-zero [4]). 

This has also importance for a system in a pure state, but if one is interested 
into the entanglement between two non-complementary parts A\ and A 2 . Indeed, 
generically the union A\ U A 2 is in a mixed state with density matrix pAiuA 2 = Tr^p 
with B the complement of A\ U A 2 . As a matter of fact, also the converse is true: any 
mixed state can be obtained by tracing out some degrees of freedom from a properly 
defined larger system, a procedure called purification (see e.g. [5]). 

A proper definition of the bipartite entanglement for a general mixed state 
(or equivalcntly the tripartite entanglement in a pure state) has been longly a 
problem because most of the proposed measures rely on algorithms rather than 
explicit expressions and are therefore hard to evaluate analytically (see e.g. Rcfs. 
[1, 4, 6, 7]). However a computable measurement of entanglement, called negativity, 
has been introduced in a seminal work by Vidal and Werner [8]. They also showed 
how the negativity provides bounds on a few operationally well-defined measures of 
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entanglement (such as distillable entanglement or teleportation fidelity). The precise 
meaning of the negativity in quantum information has been indeed established in Rcf. 
[8], but the main reason for its success is a practical one: the negativity is obtained 
from any mixed state by computing a partial transposition and then diagonalizing a 
matrix, while we still do not know generically how to compute any other mixed-state 
entanglement measures. 

Following Ref. [8], the negativity is defined as follows. Let us consider a density 
matrix p corresponding to a given mixed or pure state acting on a Hilbert space %. 
Let us consider the bipartition H = Hi ® %2 and let us denote by |ej^) and |ej 2 ') two 
arbitrary bases in the Hilbert spaces of each part. The partial transpose (e.g. with 
respect to the second space) of p is defined as 

(er ) ef ) |p^|e«ep)) = ( e « e p)|p|el 1 )ef), (4) 

and then the logarithmic negativity as 

£ = ln||p T2 || =lnTr|p T2 |, (5) 

where the trace norm ||p T2 || is the sum of the absolute values of the eigenvalues of 
p T2 . In Ref. [8] another measure, termed simply negativity, has been also introduced 

llf T2 ll - 1 

N = llP I , (6) 

which is trivially related to £ as N = (e £ — l)/2. However, £ is additive, while N is 
not and for this reason we will concentrate in the following on £ from which Af can 
be trivially derived. When the two parts are two microscopic degrees of freedom (e.g. 
spins), the negativity is equivalent to other commonly used entanglement estimators 
such as the concurrence [1, 9]. However the definition of the negativity is more 
appealing because it is basis independent and so calculable by quantum field theory 
(QFT) which naturally unveils universal features, in particular close to a quantum 
critical point. For ID critical theories, that at low energy are also Lorentz invariant, 
the powerful tools of conformal field theory (CFT) can be applied. 

The fact that the negativity is computable made it a remarkable tool to study the 
tripartite entanglement content of many body quantum systems both in their ground- 
state [10, 11, 12, 13, 14, 15] or out of equilibrium [16, 17]. Some studies for the 
bipartite entanglement at finite temperature have been presented as well [18, 19, 20]. 
However, only in a recent short communication [21], we carried out a more systematic 
and generic approach to negativity based on QFT (and in particular CFT). In the 
following we give detailed derivations of all results announced in Ref. [21] and report 
some other new findings. 



1.1. Setup and notations. 

In this manuscript we only consider the entanglement, measured by the negativity, in 
the ground state of one-dimensional systems (ID), although some derivations have a 
more general validity. We focus on the tripartitions depicted in Fig. 1, i.e. the entire 
system is divided in three parts A\, A 2 and B and we consider the entanglement 
between A\ and A 2 , whose union is generically denoted with A. As usual, A's reduced 
density matrix is pa = Tig p. In order to lighten the notation, we will denote as pi 
and pi the reduced density matrix corresponding to A\ and A 2 respectively, i.e. 

Pi = pAi = Tr^ 2 (p A ) = Tr BuA2 (p) , (7) 

92 = PA 2 = Tr Al (p A ) = Tr BuAl (p) ■ (8) 
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Figure 1. The three main configurations of ID systems we consider. Top: the 
entanglement between two disjoint intervals A\ and A2 embedded in the ground- 
state of a larger system formed by the union of A\, A2 and the remainder B. The 
whole system can be either finite or infinite. Middle: The entanglement between 
two adjacent intervals in a larger system. Bottom: The entanglement between 
two adjacent intervals which form the full system (B — ¥ 0). In all cases we denote 
A = A x U A 2 . 

Also the partial transpose with respect to A\ and A2 degrees of freedom will simply 
be denoted with the superscripts T\ and T2 respectively, i.e. 

Pa =Pa > and Pa = Pa ■ ( 9 ) 
2. A replica approach for the negativity 

2.1. Replicas and entanglement entropy. 

One of the most successful approaches to calculate the entanglement entropy for a 
bipartite system is based on the replica trick [22] which proceeds as follows. One first 
calculates the traces of integer powers of the reduced density matrix 

Trp?=^Cr, (10) 

i 

where Q are the eigenvalues of the reduced density matrix p2- Then, if one is able to 
analytically continue this expression to general complex n, the entanglement entropy 
is given by 

Sa 2 = ~ lim lim In Tr p» . (11) 

n-*i an n-»i 1 — n 

It has been shown [22, 23] that, for integer n, Trpj is a partition function on a 

complicated Riemann surface (or equivalently the correlation function of specific twist 

fields, as we shall review later) that is analytically achievable in a quantum field theory. 

The analytical properties of Tr pV^ in the complex plane of replicas are also discussed 

in Refs. [24, 25]. 
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2.2. Replicas and negativity. 

A natural way to use a replica trick to compute the negativity would be to relate it 
to the traces of integer powers of p T2 . The trace norm of p T2 can be written in terms 
of its eigenvalues Xi as 

i \i>0 \i<0 \i<0 

where in the last equality we used the normalization ^ Aj = 1. This expression 
makes evident that the negativity measures "how much" the eigenvalues of the partial 
transpose of the density matrix are negative, a properties which is the reason of the 
name negativity. 

The traces Ti(p T2 ) n of integer powers of p T2 have a difference dependence on |Aj| 
depending on the parity of n. In the following we will always indicate an even n = 2m 
as n e while an odd one n = 2m + 1 with n a . It is understood that n e and n a refer 
always to the same thing, it is just the functional dependence of the traces which is 
different. In fact, for n even and odd, the traces of integer powers of p T2 are 

Tr(p T T= = E A " e - E M" e + E M" e . ( 13 ) 

i \i>0 \i<0 

Tr(^)- = E A r° = E i^r° - E w no ■ ( 14 ) 

i \i>o \i<a 

If now we just set n e = 1 in Eq. (13) we formally obtain Tr|p T2 | in which we are 
interested. Oppositely, setting n a = 1 in Eq. (14) gives the normalization Trp Ta = 1. 
This means that the analytic continuations from even and odd n arc different and 
the trace norm in which we are interested is obtained by considering the analytic 
continuation of the even sequence at n e — > 1, i.e. 

£ = lim lnTr^ 2 )™" . (15) 

n e — >1 

We should also mention that being (p T2 ) T = p Tl , it trivially holds Ti(p T2 ) n = Tr(p Tl ) n 
for any n, even not integer. 

At first, this approach can seem rather unnatural because what we are doing 
is basically propose to calculate a quantity for even numbers and at the end set this 
even number to 1 which is instead odd. However, similar replica calculations have been 
already successfully applied to other different physical problems [26, 27], showing the 
reasonableness of the approach. 

2.3. Integer powers of the partial transpose and the negativity of a pure state. 

As a first example and check of the replica trick, we consider the case of a bipartition 
of the Hilbert spaces % — Hi ® Hi of a pure state |^) with p = |vE')( v I'| for which the 
negativity is known [8]. The Schmidt decomposition, in terms of two bases |ej^) and 
\e^} in the Hilbert spaces Hi and Hi respectively, gives 

i*> = E c *i e ? )e f > ( 16 ) 

3 

where the coefficient Cj can be chosen such that Cj <E [0, 1]. We also have 
P = Y,^\efef){e^ef\ , p 2 = Tr l(p ) = E^Xefl , 

j,k 3 
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where p 2 is the reduced density matrix on "H 2 - 

The partial transpose of the density matrix is then 

^■ = E^i e ? ) 4 a) x4 1) e? ) i. 

j,k 

whose n-th power reads 

{p T2 T = E c ji Cfe i c ^ Ck * ■ ■ ■ Cfe - 



(18) 



(2)J1) W J2)J1) U (2)J1) W (2) (1), , (2) (l)w (2) (1), 

e.„ c. | . . . \e kn Cj n A e j„ e fc„ I 



= c ji c k 1 Cj 2 Ck 2 ■ ■ ■ Cj n Ck 



x l e fei e ii ) 5 hM 5 ki,j2 S 32M S k 2 ,j 3 ■ ■■$jn-i,k n h n - 1 ,j n (e jn e kn \ . 
The sequence of deltas gives a result which depends on the parity of n 

Ee-c-|e^e«)(e^e«|, 



T 2 \n 



n = n Q odd, 
n = n e even, 



from which 



n = n a odd, 



Tr(p T2 )" 



(Trp^^ 2 ) 2 , n = n e even. 



(19) 



(20) 



(21) 



Notice in particular that 

Tr(p T2 ) 2 =Tr( / 9 T2 ) = 1. (22) 

Taking the limit n e — > 1, we recover the result [8] that for a pure state the 
logarithmic negativity is the Renyi entropy of order 1/2 (cf. Eq. (2) with n = 1/2) 



£ = S 



(1/2) 



21nTr/4 /2 . 



(23) 



Taking instead the limit n a — > 1, we recover the normalization Tr(p i2 ) = Trp 2 = 1- 



3. Negativity and Quantum Field Theory 

In this section we show how to compute for a generic tripartition of a ID quantum 
field theory the integer traces of the partial transpose of the reduced density matrix 
and from these, via the replica trick, the logarithmic negativity. For concreteness we 
will only consider the tripartition of a ID system depicted in Fig. 1 with A composed 
of two parts A = A\ U A 2 = [u\,vi] U [1*2,^2] an d B the remainder, but most of the 
following ideas apply to more general cases (e.g. A\ and A% made each of several 
disjoint intervals). In order to introduce the general formalism for the negativity, we 
first review in the next subsection the path integral approach to the entanglement 
entropy [22, 23] and the use of the twist fields [22, 24]. 
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Figure 2. Top: The reduced density matrix pa of two disjoint intervals. Middle: 
Partial transpose with respect to the second interval p^ 2 . Bottom: Reversed 
partial transpose p^ 2 = Cp^C, where C reverses the order of the row and 
column indices. 



3.1. The reduced density matrix and the entanglement entropy in QFT. 

The density matrix p in a thermal state at temperature T = 1//3 may be written as a 
path integral in the imaginary time interval (0, /3) 

P({<MIK'» = Z- 1 f t)} JJ S(cf>(y, 0) - \{ 6(<f>(y, /?) - <j> m ) e~ s - , (24) 

where Z — Tre~ H / T is the partition function and Se is the euclidean action. Here 
the rows and columns of the density matrix are labelled by the values of the fields 
{4> x } at r = 0, /3. The normalization factor Z ensures Trp = 1, and is found by setting 
{fix} = {<P' X } an d integrating over these variables. In the path integral, this has the 
effect of sewing together the two edges to form a cylinder of circumference (3. 

Now let us consider the subsystem A in Fig. 1 composed of two parts A = 
Ai U A 2 — [ui,vi] U [ii2,t>2]- The reduced density matrix p A is obtained from (24) by 
sewing together only those points x which are not in A. This has the effect of leaving 
two open cuts, one for each interval (iij, Vj), along the line r = 0. In the limit of zero 
temperature, i.e. for the ground-state of the QFT, the cylinder becomes a plane as 
in Fig. 2 (top) where the two open cuts correspond to the rows and columns of pa 
and the orientation of the arrows gives the ordering of the row/column indices, e.g. 
increasing along the directions of the arrows. 

We may then compute Tr p\, for any positive integer n, by making n copies of the 
above, labelled by an integer j with 1 < j < n, and sewing them together cyclically 
along the the cuts so that <fij(x,T = 0~) = (j)j+\(x,T = + ) and <p n (x,T = 0~) = 
<pi(x, t = + ) for all x £ A. This defines the n-sheeted Riemann surface lZ n depicted 
for n — 3 in Fig. 3. Denoting with Z-]z n the partition function on this surface we have 

Trpl=^. (25) 

This expression gives the Renyi entropies in Eq. (2) for integer n and through the 
analytic continuation the entanglement entropy in Eq. (1). 

As a consequence of its locality, the partition function can be expressed as an 
object calculated on the complex plane C, where the structure of the Riemann surface 
is implemented through appropriate boundary conditions around the points with non- 
zero curvature, i.e. the partition function in a theory defined on the complex plane 
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Figure 3. The path integral representation of Trp^ gives a n-sheeted Riemann 
surface TZ n depicted here for n = 3 and A = [u±, v{\ U [v.2, "2]- 



z = x + ir should be written in terms of certain "fields" at z 
partition function (here £[0](z, z) is the lagrangian density) 



Vj and z — uj. The 



Zn n = J [d(f>]ii n ex P — J dzdz C 



l[(j)}(z,z) 



(26) 



essentially defines these fields. Following Ref. [24], it is useful to move the topology 
of the world-sheet (i.e. the space where the coordinates x, r lie) lZ n to the target 
space (i.e. the space where the fields lie). To this aim, let us consider a model formed 
by n independent copies of the original model. The partition function (26) can be 
re-written as the path integral on the complex plane 



— 



[d(j>i ■ ■ ■ d(f) n }c exp 



— / dzdz (£[(f>i](z , z) + 
Jc 



£[</>„] (z,z)) 



(27) 



we indicated the restricted path integral with conditions 



where with J c 

j j 

C Uj , Vj ■ 0i(x,O + ) = <Pi+i(x, 0~) , x € [ui,Ui] U [u 2 ,v 2 ], i = l,...,n, (28) 

where n + i = i. The lagrangian density of the multi-copy model is 

£< B >[0i, • ■ • , <l> n ](x,T) = C[4>i](x, t) + ... + C[(/) n ]{x,T), (29) 

so that the energy density is the sum of the energy densities of the n individual 
copies. Hence the expression (27) indeed defines local fields at (uj,0) and (vj,0) in 
the multi-copy model [24]. 

The local fields defined in (27) are examples of twist fields. Twist fields exist 
in a QFT whenever there is a global internal symmetry a, i.e. J dxdTC[a<j)]{x, r) = 
J dxdTC[4>](x,T). The twist fields defined by (27), called branch-point twist fields 
[24] , are associated to the two opposite cyclic permutation symmetries i > i + 1 and 
i + 1 1— > i. We can denote them simply by T n and T n 

T n = T a , cr : i h-> i + 1 mod n , (30) 

T n = T a -i , cr -1 : i 

T n can be identified with T- n and Tn 
surface along the set A made of the two disjoint intervals [ux, v±] U [vq, V2], we have 

Trpl = {Tn(u 1 )T n {v 1 )Tn(.U 2 )Tn{v2))c ■ (32) 



1 i Y i mod n . (31) 
T ,„ = 1. Thus, for the n-sheeted Riemann 
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Figure 4. Path integral representation of Tr(p^ 2 )™ = Tr(p^ 2 )" for n = 3. 

In the following the subscript C will be understood in the expectation values, if not 
differently stated. 

3.2. The partial transposition and the negativity in QFT. 

The partial transposition of the reduced density matrix pA with respect to the second 
interval A 2 corresponds to the exchange of row and column indices in A 2 . In the path 
integral representation, this is equivalent to interchange the upper and lower edges of 
the second cut in pa as depicted in the middle of Fig. 2. If we join n copies of p T A 
cyclically, we have an n-sheeted Riemann surface where row and column indices are 
reversed compared to those of a correlation function of four twist fields, as it should 
be clear from the middle of Fig. 2. This problem can be however solved very easily 
by reversing the order of the column and row indices in A 2 as in the bottom of Fig. 2, 
to obtain the reversed partial transpose p A 2 . This is related to the partial transpose 
as p C A — Cp T A C, where C reverses the order of indices either on the lower or on 
the upper cut and satisfies C 2 = 1. Clearly Tr(p^)" = Tr(^ 2 )" and so Tr(pJ)™ is 
the partition function on the n-sheeted surface obtained by joining cyclically n of the 
above p A 2 as in the bottom of Fig. 3. In this case, the order of the row and column 
indices is the right one to identify this partition function with the four-point function 
of the twist fields 

Tr(j%) n = Tr(j%) n = (TnMTnMTMTM) , (33) 

i.e. the partial transposition has the net effect to exchange two twist operators 
compared to Eq. (32). 

For n = 2, T 2 = T 2 and so 

Trp A =Tr(p^) 2 , (34) 

which also straightforwardly follows from the properties of the trace and so it is true 
for any matrix p replacing pa above. 

To replace p T A with p A 2 it has been fundamental to consider integer cyclical 
traces. The operator C enters in quantities like Tt(paPa') which is in fact the partition 
function on a non-orientable surface with the topology of a Klein bottle. This can also 
be computed using CFT methods [28]. 
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3.3. The case of two adjacent intervals. 

Eq. (33) is of general validity, but it has interesting simple and general consequences 
when specialized to the case of two adjacent intervals. This can be obtained by letting 
v\ — » u 2 so that 

Tr(pJ)" = (T n («i)7t(«2)T„(«2)) • (35) 

For arbitrary n, this expression cannot be evaluated in general, but two cases are 
easily worked out. 

For n = 2, we have 7^ = 1 and so 

Tr(pJ) 2 = (T 2 ( Ul )T 2 (v 2 )) = Trp% uA2 . (36) 

This relation has been derived here from QFT methods and it would be interesting to 
know its general validity. 

For n — 3, we have t\ — Tz and so 

MPA 2 f = (T3(U1)T3(U2)T 3 (V2)) • (37) 



3-4- The case of a single interval. 

We now specialize to a pure state by letting B — > for which Tr(p^ 2 )™ can be worked 
out in full generality considering u 2 —¥ v\ and v 2 — » u\ 

Tr( P ^r = (T^(u 2 )Tl(v 2 )). (38) 

This expression depends on the parity of n because 7^ connects the j-th sheet with 
the (j + 2)-th one. For n — n e even, the n e -sheeted Ricmann surface decouples in 
two independent (n e /2)-sheeted surfaces characterized by the parity of the sheets. 
Conversely for n = n a odd, the surface remains a n Q -sheeted Riemann surface (a part 
from a shuffling of the sheets). This is pictorially shown in Fig. 5. Thus we have 

T:(ft) n > = (<T„ c/2 MT„ c/2 (« 2 ))) 2 = (Tr^/ 2 ) 2 , 

MPaT° = (TnJu 2 )T n Jv 2 )) = TrpX , (39) 

which are the results for pure states in Eq. (21), recovered here purely from QFT and 
showing in a first simple example the effect of the parity of n on the Ricmann surfaces. 



4. Negativity and conformal invariance: general results 

In this section we explicitly calculate the negativity for those situations in CFT where 
Tr(p^ 2 )™ can be fixed by global conformal invariance and so does not depend on the 
operator content of the theory, but only on the central charge. 

Before embarking in the full negativity calculation, we should recall that Trp^, 
for the case of a single interval A = [iti,i>i], has been obtained by uniformizing the 
n-sheeted Riemann surface to the complex plane [22] . This shows that the twist fields 
transform like primary operators of dimension [22] 

*T.=^ = (40) 
where c is the central charge. This implies (£ — \m — vi\) [29, 30, 22] 

f l\ -c/6(n-l/n) c f 

Kp n A =cJ-) => S A = -ln-+c' 1 , (41) 

V (X / O CI 
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a L n = h 



Figure 5. Cut of the Riemann surface defined by {T^ (u2)T n (v2)) ■ Left: A 
typical example for n even (n = 4 in the figure) showing how the surface decouples 
in two independent parts, each with half of the sheets. Right: A typical example 
for n odd (n = 5 in the figure) showing that the net effect is just a reshuffling of 
the sheet numeration going from (1,2,3,4,5) to (1,3,5,2,4). 

where a is an UV cutoff (e.g. the lattice spacing) and c„ a non-universal constant, 
with ci = 1. After the normalization of the two-point function of twist fields on 
the plane is fixed by the constants c„, the normalizations of two-point functions in 
arbitrary geometry and all other multi-point correlation functions of the twist fields 
are fixed and universal. Furthermore, although non-universal, the constants c n satisfy 
some universal relations (see e.g. [31]). 

4-1- A single interval. 

Even for conformal invariant theories it is useful to specialize first to the very simple 
case of a reduced density matrix corresponding to a pure system, which is obtained 
by letting B — > 0. Thus, when A 2 is embedded in an infinite system, from Eqs. (39) 
and (41) we have (setting I = u 2 ~ V2) 

MPaT* = ((T ne/2 (u 2 )T ne/2 (v 2 ))) 2 = (Trp% /2 ) 2 = < /2 (-j , (42) 

and 

Tr(/^T° = (T na (u 2 )T no (v 2 )) = Trp% = c„ o (-) . (43) 

In spite of the simplicity of the above calculation, it shows one important point of the 
CFT analysis: for n — n e even, T 2 and T n have dimensions 

2 

while for n = n odd, 7^ and T„ have dimensions 

c / 1 \ 
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the same as T Uo - Thus, performing the analytic continuation from the even branch, 
we finally have 

||^||= n limTr(^)"e=c? /2 (^) C/2 => f = ^ln^+21nc 1/2 , (46) 

which, again, just tells us that for pure states the logarithmic negativity equals the 
Renyi entropy of order 1/2. Continuing instead to n a — > 1 from the odd branch we 
obtain the normalization Trp^ 2 = 1. 

Notice that although the various constants c„ oi „ e are non- universal, they are the 
same appearing in the entanglement entropies. 

4-2. Two adjacent intervals. 

Let us now consider the non-trivial configuration in which two intervals A\ and A 2 of 
length £i and £ 2 share a common boundary (let us say at the origin) as graphically 
depicted in Fig. 1. This can be obtained by letting v\ — > u 2 = in Eq. (33) and it is 
then described by the 3-point function (we set u± = —£\ and v 2 = £2) 

Tr(p^)" = (T n (-h)T 2 n (0)T n (l2)) , (47) 
whose form is determined by conformal symmetry [63] 

(T„(-4)Tt(0)T.(4)) - 4 ^%^.^, ■ m 

which has been normalized in such a way that the structure constant C r ^=2 r is 
universal (with all the lengths £j measured in units of the UV cutoff a) and can be 
determined by considering the proper limit of the four-point function, as it will be 
done in next section. 

For n — n e even, using the dimensions of the twist operators calculated above, 
we find 

TrGoJ)^ OC (f 1 f 2 )-c/6(n e /2-2/„ e )^ i +£2 )-c/6(n e /2+l/n e ) j (49) 

that in the limit n e — > 1 gives 

i^"«(^f - *=>^k^ (») 

For n — n D odd 

Tr(p^)"° oc {£ 1 £ 2 (£i+£2))-^ {n °- 1/n ° ) , (51) 

that for n a — > 1 gives again Trp^ 2 = 1 as it should. 
Notice that for n = 2, Eq. (49) is 

Tr(p T /) 2 = clC T2nT2 (£ 1+ £ 2 )- c /\ (52) 
which equals Trp\ iUA2 , as predicted by Eq. (36), if C^^^ = c^ 1 . 
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4-3. Finite systems. 



All the previous results may be generalized to the case of a finite system of length L 
with periodic boundary conditions by using a conformal mapping from the cylinder 
to the plane. The net effect of the mapping is to replace each length £j with the chord 
length (L/n) sm(n £i/L) in all above formulas. 

Thus for the case of a pure state in a finite system the generalization of Eqs. (42) 
and (43) are 

n £,-c/3(n e /2-2/n e ) 

(53) 



Tr(^)"°=c„ o ( 



C n e /2{— 



£ 



L 

TO 



ira 

— sm — 
7ra L 



l) 



sm — 
\na L 



2 In 



-c/6(n — l/n ) 



Cl/2 



(54) 
(55) 



For the case of two adjacent intervals the finite system generalizations of Eqs. (49) 
and (51) are 



Trip*)"- cx 



-sm(— )sm(— ) 



Tr(pl') n " cx 



-c/6(n e /2-2/n e ) 



L 6 . /tt£i 
77 sm 

7T 



7r^2 
L 



sm 



L . tt(£i+£ 2 ) 
— sm 

7T L 

-c/12(n -l/n ) 



leading to the logarithmic negativity 



In 



isin (^) sin (^) 



sm 



+ enst. 



(56) 



(57) 



5. Negativity for two disjoint intervals in a CFT 

For the more difficult situation of two disjoint intervals reported in Fig. 1, the 
entanglement entropies depend on the full operator content of the CFT and this makes 
much more complicated also the calculation of the negativity. 

Global conformal invariance fixes the form of the four-point correlation of twist 
fields 

c/6(n— l/n) 



{Tn{zi)Tn{Z2)Tn{Zi)T n {Zi)) = C 

with x the four point ratio 



^31^42 



221-243241-232 



T n (x,x), (58) 



X = Z f^, (59) 
231242 

and Zij = Zi — zj. Notice that 231242/(^21243241232) = 1/(221243(1 — x)). This is 
normalized in such a way that .F„(0, 0) = 1 (for x — > 0, the four-point function is the 
product of the two two-point functions normalized with c„). The function F n {x,x) 
depends explicitly on the full operator content of the theory and must be calculated 
case by case. Notice that with J : 1 {x 1 x) = 1 as a consequence of TrpA = 1- 
Now, according to Eq. (32), Trp^ corresponds in Eq. (58) to the choice 

z\ = Mi , z 2 = vi z 3 = u 2 , z A = v 2 , (60) 



Entanglement negativity in extended systems 15 
and so 

(I \l \ \ c/6(n—l/n) 

i J" 2 -" 1 " 2 -" 1 . ,) *.(*). (61) 

(wi - ui){v 2 - u 2 ){v 2 - ui)(u 2 - Vi)J 

where, since for real Uj and Vj the ratio x is real, we dropped the dependence on 
x — x. Explicitly, in this case x is 

g= fa -»=») (62 ) 

(u 2 - u\){v 2 - vi) 

that for the order u\ < v\ < u 2 < v 2 satisfies < x < 1. 

Even Tr(p^ 2 )™ is a four-point function of twist fields according to Eq. (33), but 
it corresponds in Eq. (58) to the choice 

z\ = u\ , z 2 = vi z 3 = v 2 , z 4 = u 2 , (63) 

i.e. to the exchange z 3 -H- z 4 compared to Tr//^. Thus we can write 

(/ \/ \ \ c/6(n— 1/n) 

i v 2 ~ Ul " a " t>1 r ^ 
(vi - ui)(u 2 - v 2 )(u 2 - ui)(v 2 - vi) J 

where x = is again real and now satisfies — oo < x < 0. This can be 

rewritten in terms of an ordered four point ratio 
(vi - Ui){v 2 - u 2 ) 
(u 2 - ui)(v 2 - 



(V! - ui)(v 2 - u 2 ) x 
y = 77. 77^777. = " 7 > ( 65 ) 



with < y < 1, as 



c/6(n— 1/n) 



Tr(p A 2 ) =c„ G n (y). 66) 

\(vi - ui){v 2 - u 2 )(v 2 - ui)(u 2 - Vi) J 

Equating equations (64) and (66), we can relate the two scaling functions 7F n (x) and 

G n {y) as 

Gn(y) = (l - y) c/3( "- 1/n) (^y) • (67) 

Notice that this relation does not depend on the parity of n and so, in order to have a 
non-trivial replica limit n e — > 1 for the logarithmic negativity 



£(y)= lim In g ne (y) = lim In 

n e ^l n e — >1 



(68) 



some parity effect should appear in the n dependence of the function T n (x) for x < 0. 

The first fundamental consequence of Eq. (68) is that, for conformal invariant 
systems, the negativity is a scale invariant quantity (i.e. a function only of y) because 
all the dimensional prefactors cancel in the replica limit. This has been argued already 
in the literature on the basis of numerical data [11, 12], but never proved. 

It is useful to consider the ratio 

R ( V ) = Tr (^)" y) c/3(n-l/n) ~ 1)) m 

in which all the prefactors cancel. Being Ti(pa)" = 1 and F\{y) = 1 for < y < 1, 
we have that the negativity can be obtained just by considering replica limit of this 
ratio 

£(y)=\n lim R n Jy). (70) 
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The function T n (x) which gives the negativity has been calculated for general 
integral n only for the free compactified boson [34] and for the critical Ising model 
[35]. Some other partial results are also known [36, 37, 38, 39, 40, 41, 42, 43, 44, 
45, 46, 47, 48, 49, 50, 51, 52]. However, the calculations in Refs. [34, 35] provide 
the function T n (x) only for < x < 1 and it is a non-trivial technical problem to 
extend it to the domain x < in which we are now interested. Before embarking 
in these calculations, we discuss some simple and important physical consequences of 
Eqs. (66), (67), and (68) which are highlighted by considering the limits y — > 1 and 
y — > 0, i.e. close and far intervals respectively. 



5.1. The limit of two very close intervals: y — > 1. 

If m 2 — > v\ then y — > 1~ and we should recover the previous result for adjacent 
intervals. Let us denote with e = u 2 — v\ and £j = Vj — Uj so that Eq. (66) becomes 

/ i \c/6(n-l/n) p p 

e ' M - w,ti » = 1 - J f5r f+0(e3) - (71 » 

This reduces to the result for adjacent intervals in Eqs. (48) only if 

g n (y)~g n (l-y) a ", (72) 
with a n equal to Aj-2 the dimension of i.e. 

c /n e 2 \ 

- , if n = n e even, 

6 V 2 n„ / 



(n —\ , if n = n D odd. 



(73) 



c 

12 V"" n a 

We stress that this power-law behavior for y close to 1 is valid apart from possible 
multiplicative logarithmic corrections, whose precise form should be obtained on a 
case by case basis. 

In the replica limit n e — > 1, we then have lim a n „ — — c/4, i.e. the negativity 
diverges approaching y = 1 in the universal way 

£(y) = -- A Hl-y). (74) 
For n Q — > 1, we instead have lim a„ o = as it should to obtain the normalization 
TrpJ = 1. 

Furthermore the knowledge of the function Q n (y) allows to calculate the universal 
structure constant C r ^=2 r equating Eqs. (72) and (48), and we obtain (setting e = a) 

C T 2 =g n . (75) 
5.2. The limit of two very far intervals: y — > 0. 

The limit of far intervals y — > can be worked out from the small x expansion of 
!Fn(x) carried out in full generality in Ref. [35]. This follows from a generalization of 
operator product expansion (termed short length expansion in [35]) that for the case 
of the entanglement entropies of two intervals give 

^pa = (W ; m/b) e(^J <w»)> (^) 
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where the sum is over all the sets of possible operators {<fikj} of the CFT with 
conformal dimensions (A^. , Afc^). S{kj}{ n ) ar e calculable coefficients depending on 
the correlation functions of these operators 

n 

S{fe , } (n)a(n^(e 27ry/ "))c (77) 
j=i 

All these coefficients s^ kj j(n), in the limit n — > 1 (independently of the parity of n) 
become one-point functions which always vanish in the complex plane. This implies 
the very strong consequence that the logarithmic negativity £{y) vanishes as y — > 
faster than any power. 



5.3. The function J" n (a;) for < x < 1. 

In order to calculate Tr(/rJ)™ we first need to report some known results about the 
function J-"„(x) which has been calculated for the free compactified boson and for the 
Ising model. In the following we will discuss only the case of the boson. 

For a boson compactified on a circle of radius R (a Luttinger liquid field theory) , 
the universal scaling function F n {x) for generic integral n > 1 and for < x < 1 can 
be written as [34] 

W=*», (78, 

e(o|r) 

where T is an (n — 1) x (n — 1) matrix with elements 



= — V sin 7T- /3 fc/ „ cos 



*;=i 

and 



27T — (r — s) 
n 



(79) 



F„(x) = 2 F 1 (y,l-y, 1; or), (80) 



^ F y (z) ' 

?7 is a universal critical exponent proportional to the square of the compactification 
radius R 7 while 6 is the Riemann-Siegel theta function 

0(O|r) = exp [i7r m • r • m] . (81) 

meZ"- 1 

This generalizes the result for n = 2 in Rcf. [36]. We mention that, although the 
moments of the reduced density matrix have been obtained for all integer n, the 
analytic continuation to complex n is still beyond our knowledge and so is the Von 
Neumann entanglement entropy. 

It is also worth to report the first terms in the short length expansion (76) in the 
case of compactified free boson [35] 

/ x \ a / x \ ^ a 

F n (x) = !+(—) s 2 {n)+( — ) « 4 (n) + ..., (82) 

where a = minfy, I/77] and 

^) = "fe) a E r. ( ( 83 ) 

V4n J= i [sin(7ri)] 
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Also the coefficient s 4 (n) has been explicitly calculated [35]. Both s 2 (n) and S4(n) 
can be identified with the contributions of the short-length expansion coming from 
the two- and four-point functions of the most relevant operators in the theory. In the 
limit n4l, independently of the parity of n, we have S2(l) = s 4(l) = as it should, 
given the general considerations above. 

5.4- The function Q n {y) and the negativity for a non-compactified boson. 

Before considering the case of the compactified free boson it is worth to study the 
simpler and physical limit when the compactification radius diverges. In this case, the 
limit of Eq. (78) for r\ >• 1 has been worked out in Ref. [34] and it reads for < x < 1 

„(n-l)/2 

^n = °° (*) = 172 ' ( 84 ) 

[rife=l F k/n (x)F k/n (l - x)] 

Despite the relative (compared to Eq. (78)) simplicity of this form, many physical 
ingredients for the calculation of the negativity are already present here. 
From Eq. (67) we can write G n (y) as 

(n-l)/2 

QT°°{y) = (i - y ) ( "- 1/n)/3 v - 



1/2 : 



and the ratio with J"„(y) 

Rr°°(y) = (l-y) {n - 1/n)/3 



U , k=l F k/n{y)F k/n (l-y) 
n^i 1 Re(F fe/n (^ T )i? fc/ „( T ^)) 



1/2 



(85) 



(86) 



in which the r\ dependence drops out. The arguments of the hypergeometric functions 
in the denominator are not in the interval [0, 1] and they should be handled with care 
to select the proper analytic continuation valid for x € [0, 1]. This can be achieved by 
by using that for x € (0, 1) and k/n ^ 1/2, it holds (see e.g. [53] vol. I, section 2.9, 
Eqs. (34) and (3)) 

T(l - 2k /n) „ (k k 2k 

T(l-k/n) 2 ' 



e M / n (l-x) k / n 2 F 1 ( , k ; 

\n n n 



Fk / n (ihc) 

F k/n (^ T ) = (1 - x) k,n 2F 1 (k/n, k/n; 1; x). 



;1 



^ + (k -> n - k), 

(87) 
1/2 



We stress that these expressions are symmetric under k <-> n—k. Instead for k/n 
since we have 

F 1/2 (z) = -K(z), 



(88) 



where K(z) is the complete integral of the first kind, we can use the simpler relations 
1 



/2 



1 



= 2 -VT 

7T 



■[K(l 



iK(x)], 



(89) 
(90) 



F 1 J-^-) = -VT^K(x), 

' \X — 1/ 7T 

^^1/2(^1)^1/2(1^)] = (l)\l-x)K(x)K(l-x). (91) 

From Eq. (91) it is clear that the ratio R^^iy) equals 1 identically, as it should 
according to Eq. (34). 
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For n > 2, the ratio (86) displays a key mathematical difference between odd 
and even n. Indeed, each term in the product in both numerator and denominator is 
invariant under k <-> n — k (even if this is not evident from the formulas). Thus, for 
odd n the products in the numerator and in the denominator are both squares of the 
products for k up to (n— 1)/2. For even n this is not true because of the presence of the 
term k/n — 1/2. Nevertheless a simplification occurs between the terms k/n = 1/2 in 
the numerator and in the denominator: 

F 1/2 (y)F 1/2 (l - y) 1 
Re[F 1/2 (^)F 1/2 (^- y )] 

which is equivalent to R 2 (y) = 1. Thus we have 



(92) 



n t 'i^/n(y)Jli/ n (i-y) 

n£R«[^/n(i^T)^/n(T^ 



odd n, 



(93) 



(l-y)s("-s)-3. 



r£-l 



even n, 



ni=l R*[J , fc/n(i^l) jr fc/n(T^)] 

Despite of the relative simplicity of this formula for R n = °° (y) , we did not manage 
to work out the full analytic continuation to n e — > 1. However, it is extremely 
instructive to look at the limit y — > 1 when we can perform the analytic continuation 
explicitly and understand how the "mode" k/n — 1/2 is responsible of the main 
qualitative differences between n e even and n odd. 



5.4-1- The analytic continuation in the limit of close intervals. In order to calculate 
the analytic continuation to obtain the negativity for y — > 1~ , instead of considering 
the ratio R^°°{y), it is easier to study the function Q n {y) in Eq. (85) that we rewrite 

as 

9 n (y) = v [n - 1),2 (i-y) imn - 1/n) D n {y)-\ 

with 



D n {y) 



n— 1 ^ 

,fc=i ,y y 



1/2 



(94) 
(95) 



Each term in the product is invariant for n <-» n — k, thus we can rewrite it as 

"(n -l)/2 

n Re 



DnM 



k=l 



DnM = -(l - y) 1/2 VK(y)K(l-y) 



rie/2-1 



fe=l 



II Re[F fc/ne (-^ T )F fc/ „ c ( T ^-) 



2/ 



2/ 



(96) 



,(97) 



where for n = n e even we isolated the term k/n = 1/2 using Eq. (91). 

In order to calculate the asymptotic behavior as y — > 1~, we rewrite each factor 
in the two products above as 

MFk/niy^Wk/ni^Ty)} = (i-j/) fe/ " 2 Fi(fc/«,fc/«;i;y) 

' -r(l-2fc/n) fc/ * fc 2fc 

s n' n' n 



(98) 
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Then, for the asymptotic behavior of D n (y) close to y = 1 we need 

lim 2Fi(k/n, k/n; 1; x) — r — r th—, — < ^ , (99) 

-i- 2 sin(27rf)r(l-|)2r(^)' n 2' 

lim if(a;) = -^ln(l-a;)+21n2 + o(l), (100) 

X(0) = |, 2 Fi(o,6;c;0) = l. (101) 

For (1 - y) <C 1, in Eq. (98) the term (1 — y) k / n dominates compared to (1 - y) 1 ^ k / n 
(we recall k/n < 1/2) which can be neglected. 

Thus for odd n, the expansion of (96) as y — > 1 is 

K>-l)/2 ^p/j _ 2fcx 

D„.( S) =(i n 2sia( ^ )r(1 _"i ma) ^- »«) 

where we used Eq. (98) and JJ 1 ^ 1 ^ 2 x 2k ^ n = x^™ -1 /™)/ 4 . After some algebra, we 
have lim D no (y) = 1, as it should. Plugging the value of D no (y) into Eq. (94), we 

recover the general expansion for y close to 1 in Eqs. (72) and (73). 

For even n, using again Eq. (98) and rife^i" 1 x 2k l n = a;"/ 4-1 / 2 , we have 

D ne (y) = (1 - y) n '/ 4 ^K{y)P nc +..., (103) 

where we defined the constant 

P =1 /Tn " r "-¥> (10 4) 

and we left i^(y) instead of considering its expansion in Eq. (100), for simplicity in 
the notations. Plugging D ne (y) in Eq. (94), we recover the general expansion for y 
close to 1 of Qn e {y) oc (1 — y) c / 6 (™<=/ 2 - 2 /"e) U p £ a logarithmic correction. 
For n e — > 1, we then have 

A(y)-(i-?y) 1/4 v / ^(y)^i + o(i), (105) 

which gives the desired expansion of the logarithmic negativity for y — > 1~ 

£(y)= -ln£) 1 (y) = -^ln(l-y)-^lnX( 2 /)-lnP 1 + (l) (106) 



\ ln(l - y) - * In ( - \ ln(l - y)) - ln^ + o(l). 



(Notice a typo in the published version of Ref. [21] where we wrote the expansion of 
e £(v) instead of £ (y)). A part from a subleading double-logarithmic correction, this 
agrees with the general expansion given in Eq. (74) for the negativity of two intervals 
when they get closer. Here we can also calculate analytically the constant Pi. 
From the definition of P n we have 

— «/2-l ^ rn 2fe\ n/2-1 

2 "~ 11 2sin(^)r(l- *) 4 r(^) 11 2sin(7r^)r(l- *) 4 ' 1 ' 

where we used YYk=i 1 [r(l — 2k/n)/T(2k/n)] = 1 because in the product the same 
terms appear in opposite order in the numerator and in the denominator. Part of the 
remaining product is trivial and indeed we have 

^"fi 1 TT _ n (n-3)/2 1 2 

tt 11 2sin(7rfc/n) I 2 J n l/2 2 (l-n)/2 ^ uo > 

k=l 
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The second piece is instead more complicated. Let us consider the product 

n/2-l n/2-l 

Pn= ]J r(i - fc/«) ^ lnp„ = hr(l-ft/n). (109) 
fc=l fe=l 
Now we can employ the following integral representation of the logarithm of the T 
function 

r°° dtp- 1 / p-c*- 1 )* - i \ 
mrw-jj ^(l^—t + .-i), (no, 



to rewrite 



which, inverting the order of sum and integral, can be easily summed up as 

1 /e*/ 2 -l n\ n-2l 



r°° dtp-* 

mp„ 



o 



f 



1 - e-* Ve*/" - 1 2 



whose analytic continuation at n = 1 is 



r°° die"' 

Pi = exp ; , 



1 /e*/ 2 -l 1\ 1 

+ T, 



1 - e-* V e* - 1 2 



A 3/2 



2 l/24 e l/8 7r l/4 ^ 



(112) 



(113) 



where A = cxp(l/12 — £'(— 1)) = 1.2824... is Glaishers constant, related to the 
Riemann zeta function Q{z). Putting together Eqs. (107), (108), and (113) we finally 
have 

_2_ _ 2 7 / 6 e 1 / 2 

TTpl ~ A 6 



p i = — = 46 = 0.832056.... (114) 



5.5. Thp function G n (y) and the negativity for the compactified boson. 

As discussed in Ref. [34], for a free complex boson whose real and imaginary 
parts are compactified on a circle of radius R when encircling a branch point, 
we have the additional freedom of winding around the circle in the target space, 
i.e. (f>j(p 27Tt z 7 p~ 27TI z) — <fij-i(z,z) + 27ri?(mj ) i + irrij^), where m,j,\,mj,i € Z, 
if the branch point is in the origin. The integer numbers my can be organized 
into vectors such as mi = (mi .i, m 2j i . . . m n i) or even in 2n-dimcnsional ones 
m = (mi.i, m2.i ■ • ■ 77i„,i, 777i 2, 7ti 2! 2 • ■ • 77i„ i2 )- In all this section we will denote vectors 
with bold symbols to be easily distinguished from scalars. 

The function T n {x,x) appearing in Eq. (58) for a real scalar field is the square 
root of the same quantity for a complex field because the latter is the sum of two 
independent real scalar fields. Then, following [34] and using some results in [54, 55], 
we have that the square of T n (x, x) for a real scalar field and for any four-point ratio 
x, x £ C is given by 

n const n 

[Tn{X ^ ]2 = S F^)h,^-^h,n^)F k , n {l- X ) ^ n e_S!/ -. ( 115 ) 

where the sum runs over the possible integer components of m^, m^ 2 ' e (Z + iL) n = 
Z 2n and where the constant does not depend on x and it will be fixed later on. The 
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dependence on the compac 
first calculated in Ref. [54]: 



dependence on the compactification radius is encoded in the "classical action" Sf. 



22 

:1 

k/r, 



l _ 2gn sin(7rfc/n 



Ffc/r. 



16 



&k/r. 



(667 + 667) 



161 s 



(116) 



Pk/n ' Pk/n^ "* ' ' Pk/r, 

where g is the Lagrangian coupling C[(f>j] — g\V (j)j\ 2 / (in) and 7 = —i e -^ k / n an d £ p 
(p = 1, 2) are 



t p = 2TTRj2(9k/n) l (m i l P >+im ( l l 



7fe/r 



^2-Kik/n 



(117) 



z=i 



The dependence on the harmonic ratio x in Eq. (59) is instead encoded in the quantity 
r k / n also derived in Ref. [54]: 

F k/n (l-x) 



Tk/n 



a k/n + ifik/n ■ 



(118) 



Fk/n{x) 

being a k / n and f3 k / n the real and the imaginary part of r k / n respectively. Before 
embarking in the main calculation, we stress that the crucial difference compared to 
the case of real x with the constraint < x < 1, considered in Ref. [34], is the presence 
of the linear term in 7 and 7 in Eq. (116) which was vanishing because a k / n — for 
< x < 1. 

The products involving £ p in Eq. (116) can be written as 

66 



y^(C*fc/n)r 



(2-kR) 2 ~ Z^^ K / n ' rs ^ m rA m sS + m r,2 m s.i\ ~ { b k/n)rs [m r2 m sl - m rl m s2 \ (119) 

^ ^ r,s — 1 



»™(<?)l 



»™(<z) 



+ 



* [(&/„)„ [mgmg + mgmg] + (C- fc/n )„[mgmg - mgmg] 



where we introduced the n x n real matrices Cfc/„ and Sfc/n whose elements are 



(Cfc/n) 



^2?r^(r - s)) , (Sk/ n )rs = sin ( 27r ^( r - s )) 



cos 



(120) 



with r, s = 1, ...,n. Notice that C k / n is symmetric and Sfc/n is antisymmetric. 
Plugging (119) into (116) and summing over k, as required in (115), we obtain 



^ S t/n 



8gR tt 



2^3 



k=i 



r,s=l k=l 



Pk/n 



(^){Sk/n) rs (m^lmfl +m i r ] ) 2 m < s 2 , 

^(c k/n ) rs ( m % m % + m % m ty 



+ (C-fc/n)„ K>$ +m™m™) \ , (121) 

where many terms of S k ^ n vanish in the sum over k because they are odd under k o 

n — k. Introducing the 2n vectors of integers rhi = (m^, m^) and rh 2 = (m^ 1 ', m^ 2 ') 
made respectively by the real and the imaginary parts of and rx^ 2 \ it is possible 
to rewrite Eq. (121) as 



n 

Sf/ n = inn (rh\ ■ G ■ riii + mji • G ■ m 2 ^ , mi , rh 2 



7 2n 



(122) 



k=l 
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where n oc R 2 . The matrix G is symmetric, purely imaginary and it can be written in 
terms ofnxn block matrices as 



with 

A = 1 J2 sin(7rfc/n) C k/n = J2 sin(7rfc/n) E k/n , (124) 

n fc=i P k / n fe=i ^ fe /" 

1 ™ 1 ™ 1 
^ = ~ sin(7rfc/n) C k/n = ^ sin(7rfc/n) E k/n , (125) 

1 * ^ OL k / n r - 1 

W = — — y j — sin(7rfc/n) sin(nk/n) C k / n + cos(nk/n)S k / n (126) 

n fc=1 /4/n J 

= — fc ^" sin(7rfc/n) [sin(7rfc/n) — i cos(7rfc/n)] E k i n , 

where we introduced the complex matrices E k / n for k = 1, . . . , n whose elements are 

(E k/n ) rs EE e *ri(fc/")(r-.) /n . (127) 

These are hermitian matrices satisfying E p / n ■ E q / n = 8 p ^ q E p / ni therefore they are 
projectors, a property which guarantees that A, B and W commute. Since A and B 
are symmetric, G is symmetric by construction. The eigenvalues of A, B and W are 
respectively 

a k = sin^*), 6fc = J- sin(7^), w k = sin^)^ 7 ", (128) 

Pk/n n fi k/n n (i k / n n 

where k = 1, . . . , n. Most of these eigenvalues are degenerate because of the symmetry 
k o n — k. Since A, B and W commute, the characteristic polynomial of G is 

P(X) = det [(I - \l n )(B - AI„) - W • (129) 

(n \ n 

^2 [( a k - A)(6 fe - A) - \w k \ 2 E k/n I = Y[ [( fl fe - A )( fo fe - A) - |w fe | 
fc=i 1 / fe=i 

which is obtained from the fact that the matrices E k / n are projectors. Then the 
eigenvalues of G are 

A^r = \ (a k + b k ± v / K-6fc) 2 +4K| 2 ) (130) 

= (\r k/ n\ 2 + 1 ± ^(\r k/n \ 2 + ir-^ k/n ) ■ (131) 

They are all real and strictly positive for k e {1, . . . , n — 1}, while A* = 0. We also 
have A^ = A^_ fe for the invariance k o n — k of (128). Thus, is always a doubly 
degenerate eigenvalue. As for the degeneracy of the remaining positive eigenvalues, if 
n — 1 is even then they are all doubly degenerate as well; instead if n — 1 is odd then 
they are all doubly degenerate except for A^ 2 which are non degenerate. 
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From the common eigenvectors of A, B and W, we can construct the 2n x 2n 
matrix U diagonalizing G 



( 



G = U ■ 



V 



\ 



A 



n-l 



• u- 1 . 



(132) 



/ 



Since G has zero eigenvalues, the sum over the vectors in Eq. (115) is divergent. 
Thus, as done in [34], we introduce a regularized G e by replacing the two zeros on the 
diagonal of the diagonal matrix in (132) with a small cutoff e > 0. This allows us to 
write the regularized sum in (115) in terms of a Riemann Siegel theta function 

i-KT]{m\ G e -mi +m* ■ G e ■ m 2 ) = 6(0\r)G e ) 2 , (133) 

mi ,rri2 £ 

where the regularized G e can be written as 

e 



E 



exp 



G f 



(134) 



n \ u„ i„ 

where /„ and 0„ are nxn matrices with all equal elements given by 1 and respectively. 
Now we can isolate the divergent factor by observing that the limit 

lim [neO(0\r)G e )] = 9(0\r}G) , (135) 

gives a finite result. Since the divergent part is independent of x, it can be adsorbed 
into a normalization constant. The 2(n— 1) x 2(n — 1) symmetric matrix G introduced 
in the r.h.s. of (135) is 



G = 2i 



A W 
W t B 



where 



A = E 



n-l i |2 
T k/n\ 



Pk/n 



sin(nk/n) C, 



k/r, 



k 

n-l 



E a 7 " sin(7r/c/n) £ fe/n , 
kTi 



B = - E sin(7rfc/n) C fe /„ = V] — — sin(7r/c/n) E 1 *, /„ 

« Pfc/n. . ' Pk/n 



n-l 



Pk/r, 



W = k ^ n sm(irk/n) \ sin(TT k/n)Ci./ n + cos(7r/c/n)5fc/ n l 



k 

n-l 



T[ Pk/n 
V-k/n 



(136) 

(137) 
(138) 
(139) 



— E ~B — ~ sin(7rfc/n)[sin(7rA;/n) — i cos(7r/c/n)] E, 
k=i P k ' n 



k/r, 



These matrices are obtained from the corresponding tilded ones given in Eqs. (124), 
(125) and (126) by removing the last column and the last row, namely 



(Ck/n) rs = (Ck/n) rs , (Sk/n) rs = (Sk/n) rs , ( E k/n) rs = (-Efc/n), 



(140) 
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with r,s = 1, . . . , n— 1. Absorbing the divergent part into the normalization constant, 
we can finally write (115) as 

[T n {x,x)] 2 = const = 6(0|7?G)2 _ . (141) 

U n k Zl[F k/n (x)F k/n (l-x)+F k/n (x)F k/n (l-x)] 

As first check, for real x we have a k / n = for every k and therefore the off diagonal 
blocks W in G are zero. In this case the Riemann theta function <d(0\nG) factorizes 
into a product of two Riemann theta functions and (141) reproduces the result of [34], 
as it should. 

As for the normalization constant in (141), since it is independent of a;, it is given 
by the one computed in [34], namely (2n) n ~ 1 . Thus, for the real compactified free 
boson, we find that 

\/UkZlMFk/ n (x)F k/n (l-x)) U n k Zl\F k/n (x)yim(r k/n (x)y 

where we used that Re(F k / n (x)F k / n (l — x)) = \F k / n (x)\ 2 lm(T k / n (x)). In the limit of 
non compactified boson i] — » oo the Riemann-Siegel theta function does not contribute 
because <d(0\nG) — > 1, recovering the results for the non-compactified boson. Notice 
also that for f] = 1, while for x real and x G [0, 1] we have F n (x,x) = 1 [34], this is 
not true for general complex x and for x real and negative. 

From Eqs. (67) and (142) we can write G n (y) for < y < 1 as 



n-l 



gjy) = (1 - y )("-V«)/3 II e(o|^G(^-)y(143) 

Unfortunately, we are not able to compute the analytic continuation of C?„ c (y) to 
n e — > 1 which gives the logarithmic negativity for the free compactified boson. This 
reflects a similar problem for the standard Renyi entropies which are also known only 
for integer n > 2. 

5.5.1. Invariance of T n (x,x) for x <H> 1 — x and for r\ o 1/rj. For real x, it turned 
out that F n {x) was invariant under x -H- 1 — x. In order to show this invariance also 
for complex x, we focus on Eq. (141) in which the denominator is clearly invariant. 
As for its numerator ©(O^G), it is also invariant because when x o 1 — x we have 

1 Ctk/n Pk/n M ... 

Tk/n ^-^i- n ~ ak,n ^-KiJ^ Pk, ^KiJ^ (144) 

leading to 

A -> B, W -)• -W. (145) 

This change does not modify the sum defining Q(0\nG). Indeed, it can be reabsorbed 
by a redefinition of the sign of the first half (or, equivalently, of the second half) of 
the vector of integers over which we are summing. 

The function T n {x,x) in (142) is also invariant under r\ -H- l/rj. Indeed, by 
employing the Poisson summation formula, we can rewrite (142) in a form where this 
invariance is manifest, namely 

F n (x,x)= (146) 

Ilfc=l \F k /n{x)\ 
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T = 



where the 2(n — l)x2(n — 1) symmetric matrix T is 

irjl 1Z 
K iX/r) 

and the (n — 1) x (n — 1) symmetric matrices X and 1Z are 

~ n—1 n—1 



(147) 



X = — ^2 Pk/n sin(7rfc/n) C fe/ „ = 2 /3 k/n sin(7r/c/n) £ fe/ „ , (148) 



k=l 

n-l 



fe=l 
n-l 



ft = - ^ a k/n sm(nk/n) C k/n = 2 ^ a k/n sin(7rfc/n) E k/n . (149) 

72 k=l k=l 

Writing the quadratic form in the exponent of the generic term of the sum defining 
6(0|T), it is straightforward to see that (146) is invariant under i)H 1/n. The matrices 
(148) and (149) are respectively the imaginary and the real part of the matrix 

n—1 n—1 



T = 1Z + iX= - T k / n sin(7r/c/n) C k / n = 2 r k/n sin(7rfc/n) E k/n . 



k=i 



k=i 



We can write 9(0|T) in (146) also as 

6(0|T) = ^ exp [ Mr(p* • t • p - p* • f • p) 



(150) 



(151) 



p,p 



where 



(P,p) 



n m-y/277 n 
2r? + 2 ' 



~2 



n,m e J/ 1 - 1 . (152) 



Comparing with the results in the literature [56], from (151) we conclude that r is 
the period matrix of the genus g = n — 1 Riemann surface we are dealing with. For 
x G (0, 1) we have that a k / n — and therefore all the elements of the matrix 1Z vanish 
identically. In this case 0(O|T) = Q(0\irjX)Q(0\iX /rj) and the ?j f> \jr\ invariant 
expression found in [34] is recovered. 



5.5.2. The special case n = 2 for the compactified boson. It is instructive to analyze 
the details of the case n = 2 with complex x. Now the Riemann surface is a torus and 
the period matrix reduces to a number: its modulus ti/ 2 (x). First we recall that 



F 1/2 (x) = -K(x) = 6 3 (t 1/2 ) 2 , 
where K(x) is the complete elliptic integral of the first kind and 



n/2(x) 



.K(l-x) 
1 K(x) 



h(Tl/2 
M r l/2 



(153) 



(154) 



For n = 2 we have that T 2 = T 2 - From Eqs. (58), (146), and (153) we find that the 
four point function of the twist fields T 2 for the compactified real boson reads 

1/4 9(0|T) 



(r 2 (zi)r 2 (z 2 )r 2 (z3)r 2 (04)) 



^31^4 2 



^21^32^41^43 



mri/2) 2 



Z2lZ 3 lZ4lZ 32 Z i2 Z4 3 [x(l - x) 



1/12 



e(o|r) 



(155) 
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where T is the 2x2 matrix given by (147) for n = 2, namely 
ivPi/2 a l/2 

From (154) and the identity Q\ — Q\ — 9 2 one can find that 

The transformation properties of the Jacobi theta functions lead also to the following 
identities (here r = T1/2) 

1 



(156) 



(157) 



.t(-1/t) = 1 - x(r), 



z(t) ; 



x(t + 1) = 



x(r) - 



(158) 



Plugging (157) into (155), we find 
(72(21)72(22)7^(2:3)72(24)) 



221231241232242243 
1 



221231241Z32242243 



1/12 0(O|T) 

|e 2 (T 1/ 2)03(ri/2)e 4 (Ti / 2)| 2 / 3 

0(O|T) 



1/12 



|2 ' 



(159) 

2 1 / 3 h(ri /2 )r 

where in the last step we used the identity 2r/ 3 = 62O384 and 77(2) is the Dedekind eta 
function, which should not be confused with the compactification parameter. Now we 
observe that Z 2 {ti/ 2 ) = © (0 1 T 1 ) / 1 77 (tx / 2 ) 1 2 is the partition function of the compactificd 
boson on the torus (see e.g. Eq. (10.62) of [63] with ^hcrc = ^there/^)- 

Let us consider the exchanges 21 -R- 23, 2 2 <R- z 3 and z 3 -R- z 4 separately. From 
(59) it is easy to see that they correspond to the following involutions: x o 1 — x, 
x «-> 1/x and x x/(x — 1) respectively. From (158), one notices that they are also 
associated some SL(2,Z) transformations of Ti/ 2 - Since the prefactor containing the 
Zjj's in (159) is invariant under the three exchanges above and the partition function 
-Z2 (^1/2) is SL(2,Z) invariant [63], we conclude that the four point function (159) is 
invariant under z\ 23, z 2 -H- 23 and Z3 -H- Z4 separately. This allows us to find that 
the ratio defined in (69) for the compactified boson is given by 

Mv) = 1 ( 16 °) 

identically for < y < 1. However, we stress that, ultimately, this is just the fact that 

r 2 = r 2 . 



6. Systems with boundaries 

Let us consider a ID system on the semi-infinite line, say [0, 00) and a bipartition 
A 2 = [0,1] and A\ the reminder. We then have 

(1t\ -c/6(n c /2-2/n a ) 



Mpl 2 ) n ' = (rl(£)) = ci /2 {- i 

■ 



m / 9f\ -c/12(n -l/n ) 

c 21 
£ = - In — + 21nci/ 2 . 

4 a ' 



(161) 



The non-universal constants c„ are not the same as those appearing in the case of the 
infinite system c„, but are related to them via the Affleck-Ludwig boundary entropy 
[32] as discussed in Refs. [22, 33]. 
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Like for a periodic system, in the case of a finite system of length L with the 
same boundary conditions on both ends, the negativity can be obtained by a simple 
logarithmic mapping which results in the net replacement £ — > ^sm(n£/L) and we 
get 



2L . /tt£n 
— sin 
ira 



(?) 



21nc V2 . (162) 



The case of two adjacent finite intervals (or even one interval not adjacent to a 
boundary) is a two-point function in the half-plane, which has the same complexity as 
a four-point function in the full plane and then depends on the full operator content 
and not only on the central charge, as for the case of two disjoint intervals on the full 
line. However, some results can be achieved in a quite general way. 

Let us first consider the case of two adjacent intervals, the first one starting 
from the boundary, i.e. A\ = [0,^i] and A 2 = [£\ + £2] where B = [£ 1 + ^2,00) is 
the remainder. We place the spatial coordinate along the imaginary direction in the 
complex plane, while the imaginary time is on the real direction, i.e. z — t + ix, thus 
the system lies on the upper half plane (UHP). The traces of integer powers of the 
partial transpose reduced density matrix, correspond to the UHP two-point function 

Tr(pJ) n -(r'(zi)r„(z 2 ))uHP, with 21=^1, =i(*i+&). (163) 

The general form of such two-point function can be obtained by images method and 
it is [67] 

kill T "l 2 22l T " Zl2Z 12 

where zi = Zi, Zij = Zi — Zj, y is the four point ratio built with the two points and 
their images, and J-'niv) 1S a scale invariant function of y. Calculating this function 
is complicated, but there are some general conclusions which can be drawn without 
making any calculation. Indeed, let us assume l\ = £ 2 — £■ In this case y — 1/9 is 
constant, independent of £, as the physical intuition suggests. Then we have 

=2,._ 1_ 

£ A T% £A Trl 

Using the values of Aj-2 in Eqs. (44) and (45), we have 

( e-c/Hn-l/n) nodd, 

Tr(j%) n <x{ (166) 

[ ^-c/12(n-l/n)-c/6(n/2-2/n) n evcn . 

In particular, the n dependence for n = n e even is quite peculiar, being different from 
others found before. In the limit n e — > 1, we obtain 

£ = I In £ + const, (167) 

which however is similar to what found in other circumstances. We stress that in the 
case of a finite system of length L, the result above applies if and only if one scales L 
with £ in such a way to keep fixed the four-point ratio y. In all other circumstances, 
one should use the full formula with the unknown function T\ (y) . 

As a last example we mention the case of A 2 being a single interval detached from 
the boundary and A\ the remaining part of the half- line. In such a case, pa corresponds 
to a pure system and so the general result £ — S A \j 2 ^ in Eq. (23) applies and the results 
for the Renyi entropies discussed in Rcfs. [34, 68] apply straightforwardly. 



(T n (i£)T n (i2£)) mP cx A . (165) 
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7. The harmonic chain 

In this section we present accurate numerical checks of our CFT predictions for the 
harmonic chain with periodic and Dirichlet boundary conditions. The Hamiltonian 
of the harmonic chain with N lattice sites and with nearest neighbor interaction (but 
the procedure is easily generalized to arbitrary interactions) is 

N-l 



n=0 



1 2 MiO 2 , K, 
l^Vn + -l^ln + yWn+l - Qn) 



(168) 



where periodic boundary conditions correspond to = qa (and pn = Po) while 
Dirichlet boundary conditions to go — Qn = Pa = Pn = (these are sometimes called 
fixed wall conditions) . The variables p n and q n satisfy standard commutation relations 
[q n ,Qm] = [Pn,Pm] = and [q n ,p m ] = i5 nm . 

The chain is defined by the three parameters co,K,M, however not all of them 
are essential. Indeed, we can use a canonical transformation 

(q n ,p n ) -> ((MK) 1 ^ 4 q n ,p n /(MK) 1 ^ 4 ), (169) 

to eliminate one of the three parameters. Introducing a = y/M/K then we can rewrite 
the Hamiltonian as 



N ~ 1 r i 2 1 

1 2 ftW 2 1 



(170) 



2 ujuj 2 / \ 2 

Ya Pn + ^<l n + -(q n+1 -q n ) 

In this form, it is evident that the Hamiltonian is the lattice discretization of a free- 
boson (Klein Gordon field) with lattice spacing a and mass w. Indeed in the limit 
a — > 0, N — > oo, with Na = L, we can replace 

Pn 

q n — > 4>{x), > ir(x) — <f>(x), with x = na, (171) 

a 

(satisfying [(j)(x), n(x')] = iS(x — x') with S n , n > — > aS(x — x')) and the Hamiltonian 
above reduces to the two-dimensional euclidean action 

5 = \ 1 dx I dT K^) 2 + w ^ 2 ] ' ( 172 ) 

For uj = the theory is critical and conformal with central charge the c = 1. The 
boson field 4>{x) is not compactified and so the theory corresponds to the limit n — » oo 
in the previous CFT description. 

7.1. Diagonalization and correlation functions of the harmonic chain with periodic 
boundary conditions. 

For periodic boundary conditions, the Hamiltonian (168) can be diagonalized by 
exploiting the translational invariance and introducing the Fourier transforms of the 
canonical variables, namely 

JV-l N-l 

* = 4r E «* e27 " WAr < 9 ~ fe = TXr E «« ^ lksl \ (173) 



k=0 s=0 



(and similarly for p r ) for r = 1, . . . , N. The Hamiltonian (168) then becomes 

2M 



N-l , , „, 2 



fc=0 



^ = E ( ^7PkPN-k + ^^qkQN-k ) , (174) 
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where 



AK /7rfc\ 2 

oj k = y w 2 + _ sin[— J > w, fc = 0,...,JV-l. (175) 

Notice that a;_fe = wjv-fc = uj k . As usual, we identify pN-k with and gjv-fc with 
gLfc. The minimum values assumed by the dispersion relation uVs is w = w. From 
the canonical commutation relation [(?,■, p s ] = iS rs , one finds that q k and pN-k are 
canonically conjugate (i.e. [q~k,P-k'} = i5k,k')- Now one defines the annihilation and 
creation operators as 



ak 55 v r + **) ' fll " v (*-* - ^- fe ) ' (176) 

which satisfy the algebra [a,k,ak>] = [a\.,a k ,] = and [ak,a\.,] = i5k,k'- Then the 
Hamiltonian (174) becomes 

N-i , , x 

ff=5> fc (a{a k + - • (177) 
fc=0 ^ ' 

The ground state of the harmonic chain is then the vacuum |0) of the a^, satisfying 
afc|0) = for any k. 

The two point functions (q r q s } and (p r Ps) in the vacuum are worked out writing 
(q~k,Pk) in terms of the operators and a],, finding 

^>'M N fjk™[*^\ <178 » 



JV-l 



(0|p r p s |0) = ^- Mw fe cos 



2N 

k=0 



N 

2irk(r — s) 



N 



These correlation functions can be organized in correlation matrices Q and P whose 
elements are the correlation functions, i.e. 

Q rs = (0\q r q s \0), P rs = {0\p r Ps\0), (179) 

satisfying Q • P = Ijv/4, where Ijv is the N x N identity matrix. We also have 
(0\q r p s \0) =i5 r , s /1. 

It is very important for what follows to stress that {0\q r q s \fy is not well defined 
when <jj — because cjq = 0. This is the well known problem arising from the zero 
mode, i.e. a constant translational invariant field configuration. For w > we can 
isolate the term in Eq. (178) which diverges as w — > 

1 1 ■\r^ 1 1 \2irk(r — s)i ,„„„s 

<0| M .|0> = _ + [— ^] , (180) 

k— 1 



and the remaining sum is finite also in the limit uj — > 0. For uj = 0, the zero mode 
leads to divergent expressions, thus we work at finite but small to such that ujL <C 1. 
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7.2. Diagonalization and correlation functions of the harmonic chain with Dirichlet 
boundary condition. 

For the Dirichlet boundary conditions, the diagonalization cannot be performed by 
Fourier transform because of the breaking of translational invariancc. However, we 
can simply use the Fourier sine transform, defining as in Ref. [62] 

n 2>r sin (ivr)< ^ = Vn 2>~ feSin (iv~]' ( 181 ) 

r=l v 7 k=l V 7 

(and similarly for p r ) with r = 1, . . . ,N — 1. In this basis the Hamiltonian is diagonal 



N-l , 



1 .o MCjI 



fc=i v 



where the dispersion relation reads 



/ AK ( irk ' 
w fe = Jlu 2 + — sin 2 — > lj, k = 1, . . . , N - 1. (183) 



M \2N 

The key difference with respect to the periodic case is that ojk > to ^ for 
k = 1, ... ,7V — 1. Thus, in this model all the quantities are well defined even for 
oj = 0. In terms of the creation and the annihilation operators, the Hamiltonian takes 
the form (177) but with uj k replacing uj k . 

The correlators in the ground-state of the Hamiltonian are 

N-l 

I. I 77 1,' I \ 



7rfcr\ . /irks 



(0\p rPs \0) = - }2 MCj k sin (— ) sin (— ) = V rs , (184) 
fe=i 

where r,s = 1,...,N. 

For these Dirichlet boundary conditions we can write down exact closed formulas 
for to = and in the thermodynamic limit N — > oo, when the correlators (184) become 

(0|M.|0> = g[r («(l/2 + >• + »)- *(l/2 + r - »)) , 

(0M.|0) = %g( ^±jgM _ i-a(r-W-.)i )t (185) 

where ^(2;) is the digamma function. 



7.5. The reduced density matrix and its partial transpose. 

The construction of the reduced density matrix for the ground-state of the 
harmonic chain has been detailed in a several papers with slightly different theoretical 
approaches [57, 10, 58, 59, 60, 61]. Following Refs. [58, 61], we can easily relate pa 
to the correlation matrices Q rs = (q r q s ) and P rs = (p r p s ) (derived above for periodic 
and Dirichlet boundary conditions) restricted to the part A of the system both if the 
subsystem is connected or formed by an arbitrary number of pieces. Indeed, since the 
ground-state is Gaussian in normal coordinates (provided that no normal frequency 
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vanishes) the reduced density matrix for an arbitrary subsystem A can be written in 
the Gaussian form [58, 61] 



p A oc cxp ( - < ,''','>.,)■ 



(186) 



where &j and bj are bosonic creation and annihilation operators related to the 
original operator in the subsystem by a canonical transformation. At this point it 
is straightforward to fix the value of the tj by imposing that the correlation functions 
of g's and p's in A calculated with pa equal the ground state ones (e.g. the one 
calculated above for periodic and Dirichlet boundary conditions). Denoting with P^ 
and Qa the restriction of the matrices P and Q to the subsystem A and with p'j 
(j = 1 . . . £, with £ the number of sites in A, not necessarily connected) the eigenvalues 
of the matrix P^ • Qa, i-e. 

Spectrum^ ■ Pa) = {pi ■ ■ ■ ($} , (187) 

one finds pj — \ coth ^ [58] . (Alternatively one could work with the 21 x 2£ matrix 
formed by two-point functions of position and momentum operators, called reduced 
covariance matrix. Its reduction to a diagonal form is related to the problem of finding 
the symplectic spectrum of a symmetric and positive definite matrix [10, 59].) 

From the eigenvalues pj , using the explicit form of pa above, we obtain the Renyi 
entropies 



Trp A = n 



Pj + 



and the entanglement entropy as 



Pj + 



In 



Pj 



Pj 



Pj 



In 



Pj 



(188) 



(189) 



The partial transpose of the reduced density matrix has been constructed in 
Ref. [10], where it has been shown that the partial transposition with respect to the 
subsubsytem A^ has the net effect of changing the sign of the momenta corresponding 
to the subsystem A2, leaving the partial transpose reduced density matrix p^ 2 a 
Gaussian matrix of the form (186). This means that we can simply replace the matrix 
P^4 with a matrix in which all the momenta in A2 have been changed sign, i.e. 

(190) 

where Wa 2 is the £ x £ diagonal matrix having —1 in correspondence of the sites 
belonging to A 2 and +1 otherwise, i.e. (KA 2 )rs = 5 rs (—l) SreA 2 
because R^i = — Ka 2 ■ Then the eigenvalues Vj of the product 
full spectrum of p 7 ^ . Denoting 



5T2 

■ A 



Notice that P^ 1 



IA ' I 



^ 2 determine the 



Spectrum (Qa 



we have 



Tri 



Pa 



n 

3=1 



(191) 



(192) 
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Figure 6. Entanglement in a periodic chain of total length L = 100. We consider 
the bipartition between an interval of length I and the remainder. Here z = l/L. 
Left: Tr(p T2 )" for n = 3,4, 5,6 as function of sin(7rz). The straight lines are the 
CFT predictions where the only free parameter is the overall amplitude c n (cf. 
Eqs. (53) and (54)). Right: The logarithmic negativity E vs the CFT prediction 
(55). The three curves correspond to different values of ui and show the influence 
of the zero mode. Indeed they are parallel (in log-linear scale) and the zero mode 
affects only the non-universal additive constant. 



and the trace norm 
\\p T a\ 



n 

i=i 







1 




1 






V 3 ~ 


2 




J 2 





max 



2iA 



(193) 



From this quantity we easily get the negativity Af and the logarithmic negativity £ . 
Notice that only Vj < 1/2 contribute to (193). 

Let us then summarize the practical procedure to hnd the negativity and Tr(p^ 2 )' 1 
for an arbitrary tripartion A\ U U B in the ground-state of an harmonic chain: 

• From the correlation matrices Q and P (which for periodic and Dirichlet boundary 
conditions are given in Eqs. (178) and (184) respectively), construct the reduced 
correlation matrices Qa and Pa by erasing the rows and columns corresponding 
to the part B of the system. 

• Change the signs of the momenta in the part A- 2 to construct the matrix P^ 2 as 
in Eq. (190). 

• Calculate the eigenvalues of the product Qa ■ ^a ■ 

• Use Eq. (192) to calculate Tr(p^ 2 )" and Eq. (193) to obtain the negativity from 



In the semi-analytic calculations presented in the following we will always work 
in units M — K — a — 1 and so the total length is L — Na = N. 
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sin(n:z) Lsin(:nz) 

Figure 7. Entanglement in a critical (a; = 0) harmonic chain with Dirichlct 
boundary conditions of total length L. We consider the bipartition between the 
interval Ai = [0,1] and the remainder A 2 = [£, L], Here z = IjL. Left: Tr(p T 2) n 
for n = 3,4,5,6 and for L = 200 as function of sin(-7rz). The straight lines are 
the CFT predictions where the only free parameter is the overall amplitude c n 
(cf. Eq. (161)). Right: The logarithmic negativity £ for L = 25, 50, 100, 200. 
The four sets of data collapse on the same curve when plotted against the chord 
length, as they should. 



7.4- The negativity for one interval. 

We consider here the case when A\ is a block of I consecutive sites and A 2 the 
remainder, i.e. B — > 0, as shown in the bottom of Fig. 1. In this case, the 
negativity and the traces Tr(p T2 )" coincide with the Renyi entropies for any pure 
state (cf. Eq. (21)) as proved in Sec. (2.3). It is however important to report these 
results for a twofold reason. On the one hand, knowing a priori the final result of 
the calculation provides a non-trivial check of the numerical procedure which indeed 
differs substantially when calculating Tr(p T2 )™ and Trp% since for latter there is no 
need neither of a partial transposition, nor of the correlation matrices in the part A\ 
(that involves basically the same procedure for one and two intervals). On the other 
hand, these controlled calculations give also an idea of the corrections to the scaling 
to the asymptotic CFT formulas. 

In the case of periodic boundary conditions, we cannot work directly with lo = 0, 
because of the zero mode. Thus, we consider several values of ui imposing the further 
constraint ujL <C 1 to ensure that all the data for any 1 < I < L are in the conformal 
regime. The data for Tr(p T2 )™ in the periodic harmonic chain are reported in Fig. 
6 for n = 3,4,5,6. The agreement with the conformal predictions in Eq. (54) is 
excellent. The visible deviations for small £ are the corrections to the scaling to the 
entanglement Renyi entropy discussed in Ref. [64, 65, 66] of the unusual form l~ 2 / nR 
where tir is the index of the corresponding Renyi entropy (i.e. — n for n odd and 
n R = n /% for n even). 

In Fig. 6 we report also the logarithmic negativity as function of the chord length, 
finding a perfect agreement with the conformal prediction (55). We report the data 
for three different values of u>, all satisfying ljL <C 1 and so in the conformal regime. 
It is evident that in logarithmic-linear scale the three curves are parallel confirming 
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Figure 8. A typical example of the entanglement of two adjacent intervals of 
the same length i in a periodic chain of length L. The legend refers to both plots 



and we fix u = 10~ B . Left: Tr(p^ 2 ) 3 as function of z = t/L. The data are 
consistent with the conformal prediction (full lines with only one multiplicative 
free parameter). Right: Logarithmic negativity E for the same chains as on the 
left. Again the data are perfectly consistent with the CFT predictions (full lines 
with only one additive free parameter). 



that the zero mode only affects the value of the non-universal additive constant and 
not the leading logarithmic behavior of £ with the subsystem sizes. 

We then perform the same analysis for a finite system of total length L with 
Dirichlet boundary conditions considering the interval A\ starting from the boundary 
up to £ and the remaining L — £ sites. In this case, as discussed above, there 
is no zero mode and we can work directly at u = 0. The data for Tr(p T2 )™ are 
reported in Fig. 7 for several values of n, all for L = 200. It is evident that 
increasing I the data approach the CFT predictions in Eq. (161), but in this case 
the corrections to the scaling are much larger than in the periodic case. This does not 
come unexpected, indeed in Refs. [64, 65, 68] it has been shown that in the presence of 
the boundaries the corrections to the asymptotic results are of the form (where 
again ur is the index of the corresponding Renyi entropy), i.e. they have an exponent 
which is half of the corresponding one for periodic systems. On the right of Fig. 7 
we report the logarithmic negativity as function of the chord length Lsui(-k£/ L) for 
L = 25, 50, 100, 200. All the data at different L collapse on the same curve that for 
large enough chord length are perfectly described by the CFT prediction (162). The 
corrections to the scaling are smaller than the ones for n > 2, in agreement with the 
general analysis of the Renyi entropies [64, 65, 68]. 

Finally, we checked that all the results reported in this subsection satisfy the 
relations in Eq. (21) between Tr(p T2 )™ (and £) and the Renyi entropies in the same 
state. 

7. 5. The negativity for two adjacent intervals in periodic chains. 

We now consider the case of two adjacent intervals in a periodic system of total length 
L. For simplicity we consider the two intervals to have equal length £, thus all the 
results depend on the single parameter z = i/L e [0,1/2]. In terms of z, the CFT 




Figure 9. Entanglement for two adjacent intervals of equal length i < L/2 in a 
periodic chain of length L. Left: the ratio r n (z) in Eq. (196) as function of z = l/L 
compared with the parameter free CFT prediction. Right: Subtracted negativity 
e(z) in Eq. (199) again compared with the parameter free CFT prediction. 



predictions in Eqs. (56) can be rewritten as (we fix c = 1) 

f (sin(^z))-^/ 2 - 2 / n =)/ 3 (sin(27rz))-("=/ 2 + 1 / n =)/ 6 , 

Tr(p%T <* I (194) 
[ (sin 2 (7rz) sin(27rz))-("°- 1 /«o)/i2 j 

and the proportionality constants depend on L in a known manner. For the logarithmic- 
negativity from Eq. (57), we have 

1 sin 2 (7rz) 1 
c = — in . . — r + cnst = - In tan(7rz) + cnst , (195) 
4 sm(27rz) 4 

where again the additive constant depends on L in a known manner. These predictions 
are checked against the numerical data in Fig. 8 where the only free parameters in 
each curve is fixed by a fit. The agreement is clearly excellent. 

However one can do even better checks of the theory by constructing ratios in 
which all the dependence on non-universal parameters (such as those coming from the 
zero mode) and also the dependence on L cancel. A straightforward idea would be 
to divide Tr^^ 2 )" for the value it assumes at a given fixed (., e.g. I = L/A, i.e. by 
considering the logarithm of the ratio 

r n (z)=hx TliP T { 2 2 n > ( 196 ) 
Tr(p A A2=L/i ) n 

whose parameter free CFT predictions for n even and odd are 

1/2 n„\ , 9/ . % 1 / rip 1\ 

ln(2sinVz)) - - + — J ln(sin(2^)) , (197) 

(— -n ) ln(2sin 2 (7Tz)sin(27rz)). (198) 
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Figure 10. For a periodic chain of length L, we report the ratio R n (y) defined 
in Eq. (69) as function of y for several L and for n = 3,4. The continuous lines 
are the parameter free CFT predictions to which the data converges for L — > oo. 



These ratios are shown in Fig. 9 for n = 3 and n = 4. The agreement of the numerical 
data and the CFT prediction is excellent. Some very small deviations are visible for z 
close to and 1/2, but we checked that they go to zero increasing L in a controllable 
way. 

For the logarithmic negativity, we can analogously define the subtracted quantity 

e(z) = £(£, L) - £(L/4, L) = - ln[tan(vrz)] , (199) 

and again the rhs is a parameter free CFT prediction. In Fig. 9 (left panel), this 
prediction is compared with the numerical data and the agreement is extremely good 
since no deviations are visible even close to the boundaries. 

7. 6. The negativity for two disjoint intervals in the periodic chain. 

To conclude our analysis of the periodic chain we consider now the most difficult 
situation of two disjoint intervals for which the negativity has been already considered 
numerically [12]. Here we start by considering the traces Tr(p^ 2 )" which, in the 
conformal regime ujL <C 1, should be described by Eq. (66) with Q n {y) given in (84). 
The direct numerical data, that we do not present here, agree rather well with the 
CFT where the overall constant is fixed by a fit and explicitly depends on the values 
of to and L, since the data are influenced by the zero mode. However, this is not a 
very accurate statement on the true entanglement because the function Q n {y) turns 
out to be very close to a constant (in proper units equals 1). Thus the direct data 
mainly probe the prefactor to G n {y) in Eq. (66), which however gives no contribution 
to the negativity in the replica limit n e — > 1. 

As already stated in Sec. (5), a practical way to get rid of the prefactor is to 
consider the ratio R n (y) in Eq. (69), where also the non- universal parts due to the 
zero mode cancel and we are left with a universal function of y. The CFT prediction 
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Figure 11. Finite size scaling analysis for d n (y) in Eq. (200) for n = 3 (left) 
and n = 4 (right). We report from Fig. 10 several values of y (increasing in the 
direction of the arrow, but we do not give the actual value to simplify the reading 
of the plot). The data are compatible with a leading correction to the scaling of 
the form L" 2 /™ . 



for this ratio i?^ =00 (j/) in Eq. (86) is compared to the numerical data in Fig. 10. As 
L increases, the data approach the CFT result. The differences with the asymptotic 
formula are due to the presence of unusual corrections to the scaling [64, 65] whose 
leading part is of the form L~ 2 l n . A quantitative finite size scaling analysis is reported 
in Fig. 11 for n = 3,4 showing that the difference 

d n (y)^Rn(y)-K =00 (y), (200) 

for several values of y is of the expected form L~ 2 l n . As well known (even analytically) 
for other simpler cases [64, 66] for larger n, the sublcading corrections to the scaling, 
of the form L~ 2p / n with p integer, cannot be neglected and a proper analysis requires 
the introduction of some fitting parameters. 

Finally we turn to the study of the negativity £ reported in Fig. 12 showing that 
all data collapse on a single curve, without sizable corrections. Unfortunately we do 
not have the analytic continuation of i?J' l =00 (j/) to n e — > 1 as a function of y. However 
we can study the two interesting regimes of far and close intervals corresponding to 
y — > and y 1~ respectively. For small y, the data, being very close to zero, are 
consistent with the prediction that they vanish faster than any power. In Ref. [12], on 
the basis of the numerical data, the form e~ a l has been proposed. This proposal is 
shown in logarithmic scale on the inset of Fig. 12 together with a simple exponential 
e~ b ' v '. The two possible scenarios are very difficult to be distinguished on the basis of 
the numerical calculations involving exponentially small numbers and the ambiguity 
can be resolved only by analytically continuing R r ^°°{y). 

For y — > 1, the general prediction from CFT (74) for c = 1 is 

£(y) = -\\n(l-y) + ... , (201) 

however, we have shown in Eq. (106) that for the model at hand, subleading double 
logarithmic corrections are present 

£(y) = - i]n(l -y)-\ lnK(y) -lnP 1 + o(l) , (202) 

where K (y) is the elliptic integral of the first kind and Pi — 0.832056 ... is given in 
Eq. (114). Written in this form, the expansion contains also some of the subleading 
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Figure 12. The negativity £ (y) is a universal scale invariant function 
with an essential singularity at y = 0. We report the data for L = 
50, 100, 150, 200, 250, 300 but, since they are hardly distinguishable, we do not 
give a legend box. The solid line is the expansion close to y ~ 1 in Eq. (202) 
which very surprisingly describes well the data down to y ~ 0.3. The inset shows 
the same plot in logarithmic scale showing that for small y the two possibilities 
£ ~ e~ a l y and £ ~ e~ b / are too close to be distinguished. 



corrections (in 1 — y) and it is expected to describe the data better. Indeed in Fig. 
12, this prediction is almost indistinguishable from the data all the way from y ~ 1 
(where it is an exact result) down to y ~ 0.3. We should mention that the subleading 
logarithmic correction may be responsible for the exponent 1/3 found in Ref. [12] as 
compared with our analytic result 1/4. 

Finally we would like to mention that, in a long enough chain, when each interval 
contains a finite number of lattice points, the negativity must vanish exactly for 
sufficiently large separations. This is because the reduced density matrix p& has all 
strictly positive eigenvalues. Thus, when we take the partial transpose, the change in 
the density matrix, and therefore in the eigenvalues, can be made arbitrarily small since 
(piPj) (which is the correlator that changes sign, cf. Eq. (190)) decreases like \i — j\~ 2 ■ 
Indeed this is consistent with the well-known result [1, 9] that the entanglement of 
two far away sites is exactly zero. 

7. 7. Tripartite chains with Dirichlet boundary in the origin 

Now we consider the non-trivial case of a tripartite chain on a system with boundaries 
discussed in Sec. 6, which is the semi-infinite line, with A\ = [0, £], A2 = [£, 2£] and B 
the remainder. In this case, the results for Tr(p^ 2 )™ are given in Eq. (166) which we 
report also here: 




n odd, 



n even. 



(203) 
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Figure 13. Ti(py^) n (left) and logarithmic negativity (right) for a semi-infinite 
critical (uj = 0) harmonic chain with Dirichlct boundary conditions in the origin. 
We report the results for two adjacent intervals of equal length I with the first 
one starting from the origin. The data are perfectly described by the asymptotic 
CFT predictions reported as solid lines. 



The most interesting feature is the peculiar behavior for n even, which can be explicitly 
checked by considering the semi-infinite harmonic chain with uj — 0. The results for 
the numerical evaluation of Tr(p^ 2 )™ are reported in Fig. 13 (left), showing the perfect 
agreement of both even and odd n = 3,4,5,6 where only the global amplitude has been 
fixed with a fit. The right panel of the same figure shows the corresponding logarithmic 
negativity which turns out to be described by the CFT prediction £ = l/41n^+ const. 

8. Some scaling considerations for massive theories 

It is well known that for a gapped one-dimensional model, increasing I the 
entanglement (Renyi) entropy saturates to a finite value [30]. The calculation of 
this saturation value is generically complicated because it depends on microscopical 
details of the model and so it is calculable only for few simple integrable cases (see 
e.g. [22, 69, 70]). However, some general results can be worked out when a system is 
close to a conformal quantum critical point, and in the scaling limit where the lattice 
spacing a — > and with the correlation length fixed and large. Under these hypotheses 
and when all the lengths of the various subsystems (and so the total one for a finite 
system) are much larger than £ which is itself large, the Renyi entropies are [22] 

sM=^L(l + I)log*+0(£°), (204) 
12 V n/ a 

where A is the number of boundary points between A and its complement and c is the 
central charge of the conformal field theory for £ _1 =0. Clearly, when the interval 
lengths are of the order of £, a complicated and universal (within the scaling limit) 
crossover takes place, which has been worked out only in very few QFT [22, 24, 71, 72]. 

The same difficulties for the entanglement entropy, with also some additional 
problems, make prohibitive the calculations of the negativity and Tr(p^ 2 ) n for a general 
gapped model. However, in the scaling limit and when the correlation length is itself 
large, but smaller than all the subsystems lengths, scaling considerations lead to very 
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general results. For example in the case of a bipartition of a pure state in a finite 
interval of length I ^> £ and the remainder, the scaling suggests to replace I with £ in 
Eqs. (42) and (43) for Tr^ 2 )" and so the negativity is 

£ = ^ln£ + 0(a. (205) 

This is clearly true because it is the general relation (21) between negativity and Renyi 
entropy. A direct check of this is the general result for a bisected harmonic chain (i.e. 
A2 formed by N/2 consecutive sites and A\ the remainder) with periodic boundary 
conditions and nearest neighbor interaction. In this case, an interesting exact formula 
has been found for the logarithmic negativity [10] 

1 / AK \ 

which, remarkably, is independent of the size of the chain N. In the limit £ _1 oc lo <C 1 
gives £ = (ln£)/2, in agreement with the general scaling. 

A more interesting example can be obtained by considering the case of two 
adjacent intervals. Again, scaling suggests to substitute in Eq. (50) for the conformal 
case £i t 2 with £ giving 

£= ^lne + O(e ). (207) 

This scaling in £ is checked in Fig. 14 (left) for the harmonic chain for two adjacent 
intervals of the same length £. This shows an excellent agreement for 1 <C £ <C t while 
for larger £, the crossover to the CFT prediction (50) takes place. 

On the same lines as above, the negativity for other tripartitions can be deduced 
by means of simple scaling arguments. 
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9. Conclusions 



We introduced a general field theoretical formalism to calculate the negativity and the 
logarithmic negativity, as an extension of our previous short communication [21]. This 
novel approach is based on a replica calculation of the traces of even integers n e powers 
of the partial transpose of the reduced density matrix and analytically continuing this 
to n e — > 1, i.e. the logarithmic negativity is 

£ = lim lnTr^ 2 )™" . (208) 

n e ^ 1 

Several physical situations have been explicitly worked out for a conformally invariant 
theory: 

• The case in which A\ is an interval and A 2 the remainder of an infinite, semi- 
infinite or finite (both periodic and with boundaries) system. 

• The case in which Ax and A 2 are two adjacent intervals of length l\ and l 2 
respectively. For an infinite system the negativity is 

This is simply generalized to a finite periodic system. 

• For the case in which Ax and A 2 are two disjoint intervals (always of length 
tx and £ 2 ), the negativity turns out to depend on the full operator content of 
the theory and it is a scale invariant function (i.e. a function of the harmonic 
ratio of the four points defining the two intervals). We calculated the traces of 
integer powers of the partial transposed reduced density matrix for a free bosonic 
theory, both compactified (i.e. the Luttinger liquid field theory) and in the limit 
of infinite compactification radius. However, the n dependence of these formulas 
is too complicated and we managed to calculate the analytic continuation only in 
some physical relevant limits (far and close intervals). 

All the CFT results have been accurately checked for a chain of harmonic oscillators, 
finding perfect agreement. We also proposed a few scaling relations valid in the scaling 
regime when a theory is close to a conformal quantum critical point, and the correlation 
length (inverse gap or mass) is finite but large. 

We are currently working out the generalization of the CFT approach to finite 
temperature field theories which is not as straightforward as it could naively appear. 
Finally, we mention that it is of extreme interest to check numerically all our CFT 
predictions in more complicated lattice models such as spin-chains and itinerant 
fcrmions described by the Luttinger liquid field theory, for which we worked out explicit 
predictions for Tr(p^ 2 )™ even in the case of disjoint intervals. 
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